A percolation model for slow dynamics in glass-forming materials 
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We identify a link between the glass transition and percolation of mobile regions in configuration 
space. We find that many hallmarks of glassy dynamics, for example stretched-exponential response 
functions and a diverging structural relaxation time, are consequences of the critical properties of 
mean-field percolation. Specific predictions of the percolation model include the range of possible 
stretching exponents 1/3 < (3 < 1 and the functional dependence of the structural relaxation time 
r a and exponent (3 on temperature, density, and wave number. 



As temperature is decreased near the glass transition, 
the structural relaxation time in glassy materials in- 
creases by many orders of magnitude with only subtle 
changes in static correlations [1]. In addition, structural 
correlations display an anomalous stretched-exponential 
time decay exp(— i/r Q )' 3 , where /3 is the stretching ex- 
ponent and r a the relaxation time. Understanding the 
origin of these behaviors is one of the most important 
outstanding problems in statistical physics. 

Although stretched-exponential relaxation is common 
to many glass-forming materials, the dependence of r Q 
and (3 on temperature and density is not universal. For 
molecular, colloidal and polymer glasses, where struc- 
tural relaxation is measured using density autocorrela- 
tion functions, the temperature dependence of r Q is af- 
fected by the fragility [2j. In magnetic glasses, where 
structural relaxation is measured using spin autocorrela- 
tion functions, r Q depends on details of the microscopic 
interactions [3| . In all glassy systems the stretching expo- 
nent j3 varies between 1 /3 and 1 depending on scattering 
wave vector, density and temperature, and its depen- 
dence on these variables is not universal 0]. 

How do we understand anomalous relaxation in glass- 
forming materials where, although correlation functions 
display stretched-exponential relaxation, the tempera- 
ture and density dependence of t q and /3 vary from one 
material to the next? To answer this question, it is im- 
portant to consider how the energy landscape, which is 
the energy hypersurface as a function of all configura- 
tional degrees of freedom, affects dynamics Activa- 
tion out of energy minima is rare at low temperature, 
and the system is confined to basins surrounding energy 
minima with infrequent hopping between minima 6] that 
allows the system to explore configuration space [7( . 

In this work, we focus on the connection between 
anomalously slow dynamics in glass-forming materials 
and percolation of mobile domains in configuration space. 
The decay of structural correlations over a time t is re- 
lated to the average distance that the system moves in 
configuration space during that time. Thus complete 
relaxation — decay of structural correlations to zero — 
only occurs after the system diffuses through a net- 
work of basins that percolates configuration space. For 



both hard-sphere systems with infinite energy barriers 
and systems with finite barriers, we demonstrate that 
a percolation transition in configuration space is re- 
sponsible for several hallmarks of glassy dynamics: (1) 
stretched-exponential relaxation of structural correla- 
tions and the experimentally observed range of values 
and wave-number dependence of the stretching exponent 
/?; (2) the form of the divergence of r a , and (3) a diverging 
length scale near the glass transition for hard spheres. 

A. Hard spheres: We first consider collections of 
hard spheres that interact at contact with an infinite re- 
pulsion. At moderate density, hard spheres behave as 
simple fluids. However, as density increases structural re- 
laxation becomes anomalously slow. Upon further com- 
pression, if crystallization is avoided, the system becomes 
confined to a disordered collectively jammed [8] (CJ) state 
at packing fraction <j>j. We focus on situations where, due 
to particle size polydispersity or fast quench rates, crys- 
tallization is avoided, and the systems are disordered. In 
CJ states at all single and collective particle displace- 
ments cause particle overlap; thus, due to hard-sphere 
constraints, no motion occurs at <f>j. A system of N hard 
spheres in d dimensions has ~ N\e aN CJ states 0, [l(| 
that can be represented as points in eW-dimensional con- 
figuration space, where a > is a constant. 

The transition from glass to liquid can be understood 
as the percolation of a network of bonds between points 
that represent CJ states in configuration space, where a 
bond is formed when the system can move between two 
CJ states without violating hard-sphere constraints. At 
high packing fraction there is no percolating cluster and 
relaxation is limited by the maximum cluster size. At 
low packing fraction, a percolating network of bonded C J 
states spans configuration space and the system can fully 
relax. To use the predictions of continuum percolation it 
is useful to partition configuration space into basins of at- 
traction, each associated with an individual C J state 11 1 . 
At a given packing fraction <j> we define a mobility domain 
for each CJ state as the portion of the basin of attraction 
that satisfies the hard-sphere constraints at <j>. Thus at 
= (when there are no hard-sphere constraints) the 
mobility domain and basin of a CJ state are the same, 
for cj) > the mobility domain is a subset of the basin, 
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and at 4>J the mobility domain is the CJ state. 

We formulate a percolation model by considering the 
motion of hard spheres near <f>j. At <fij the system is con- 
fined to one CJ state and no motion occurs (Fig. []Ja)). 
For (f> < 4>j there is a closed mobility domain accessible 
to each CJ state [1] (Fig. [Ub)). For smaller 4> the mo- 
bility domains of two CJ states can come into contact, 
which allows the system to transition between CJ states 
(Fig. Hljc) ) . The system will diffuse on a network of mobil- 
ity domains if many are in contact. We also expect that 
when the mobility domains become sufficiently large, a 
percolating cluster of domains is formed (Fig.QJd)). 

Structural relaxation in dense hard-sphere systems oc- 
curs via dynamical heterogeneities [1 31 ] - which implies 
that the shape of mobility domains is complex. Near CJ 
states mobility domains are roughly <i./V-dimensional and 
are quickly explored. Further from CJ states, the effec- 
tive dimension of mobility domains can be significantly 
reduced, and the time needed to explore these regions is 
large. We also expect polydispersity of mobility domain 
sizes for different CJ states. However, we assume there 
is an upper critical dimension D* of configuration space 
above which mean-field theory accurately describes crit- 
ical exponents. In this limit, complexities arising from 
the geometry and correlations of mobility domains can 
be ignored. We focus on large N where dN > D* and 
construct a mean-field theory in terms of the packing 
fraction of mobility domains II in configuration space. 
Percolation occurs at a critical value lip and is controlled 
by the mean- field percolation exponents [3] • We employ 
the mean-field assumption throughout this section, and 
outline a non-mean-field approach at the end of section B. 

Our percolation model predicts a glass for II < lip 
and liquid for II > lip. Structural relaxation of hard 
spheres can be quantified using the incoherent part 
of the intermediate scattering function (ISF) &q(t) = 
Y^j exp [iq ■ Arj(t)]/N, where Afj(t) is the displacement 
of particle j over time t and q is the scattering wave vec- 
tor [IB} . The infinite-time value of the ISF f$ = ^(oo) 
is an order parameter for the glass transition that is zero 
for a liquid and positive for a glass. For II < lip there is 
no percolating cluster of mobility domains, > 0, and 
the system is a glass since the maximum distance it can 
diffuse, which is set by the percolation correlation length 
£ cx (lip — IT) -1 / 2 , is finite. Using a cumulant expansion 
we predict log/^ oc — q 2 £ 2 oc -q 2 (U P - II) _1 for q£ <C 1. 
For II > IJp the largest cluster percolates and /§• = 0. 

Time dependence of the ISF can be determined from 
dynamics in configuration space. The model contains 
two important time scales: short times where the system 
is confined to a single mobility domain and long times 
where the system diffuses on a connected cluster of mo- 
bility domains. The short-time dynamics is characterized 
by the time To needed to transition from one mobility 
domain to another. When II > lip and a percolating 
network exists, long-time relaxation of the ISF is char- 
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FIG. 1: (Color online). Accessible regions in configuration 
space for hard spheres, (a) At <j>j only CJ states (points) are 
accessible; (b) At <j>i < cj>j motion occurs in mobility domains 
surrounding CJ states; (c) At (j>2 < <f>i transitions between CJ 
states with contacting mobility domains (shaded) occur; (d) 
At cj>p < 4>2 at least one mobility domain percolates (shaded 
yellow) and the system transitions from glass to liquid. 

acterized by anomalous diffusion on fractal percolating 
networks. The mean-square displacement (MSD) Ar 2 (t) 
on percolating networks obeys the scaling form 



Ar 2 (t)oc(t/ro) 1/3 Gfr 2 (Vro) 1/3 



(1) 



where £ cx (II — lip) x l 2 is the percolation correlation 
length and the scaling function G[z] cx z 2 for z > 1 
and G\z] = 1 for z < 1 [H, \rft. To first order in a 
cumulant expansion <E>g(t) = exp(— q 2 Ar 2 (t)/6), and the 
alpha relaxation time r a , defined via $^(r a ) = e _1 , is 



T a CX 



to q~ 



2 (U-U P ) 
T q- 6 



- 2 for g£<l, 
for q£ > 1. 



(2) 



Note that there are two contributions to t q : (1) the aver- 
age transition time between mobility domains To and (2) 
a time-scale for diffusion on percolating networks, which 
results in the factor (II — IIp)~ 2 . The ^-dependence in 
Eq. @ is consistent with experiments on dense colloidal 
suspensions [l8| and vibrated granular materials [lj|. 
Since the MSD in Eq. (fl]) crosses over from anomalous 
diffusion Ar 2 (t) cx t 1 ' 3 at small times to normal diffusion 
Ar 2 (t) oc t at large times, the percolation model predicts 
that j3 varies with time and satisfies 1/3 < /? < 1. In ex- 
periments, (3 is measured by fitting the ISF to exp(</r Q ) /3 
near r a , and thus time dependence has not been ob- 
served. Near r Q the percolation model predicts that j3 
only depends on q£ with (3 — 1/3 for q£ ^> 1 and (3=1 
for g( « 1. These limiting values for large and small q 
have been observed in experiments of hard spheres [19j . 
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To further test the predictions of the percolation model 
we express IT and tq in terms of </>. II is equal to the 
number of CJ states (N\e aN ) multiplied by the average 
volume V of each mobility domain in configuration space 
(normalized by the total volume L Nd ). V is determined 
by the total free volume Vf accessible to the system in 
real space and to first order V X ^ N — vj/(L d ) — ((f) j — 
4>)/(j)j. Thus for large N, IL = N N (l - 4>/(t>j) N . In the 
mean-field approximation, To is inversely proportional to 
the average number of mobility domains in contact with 
any single domain. Assuming a random distribution of 
hyperspherical mobility domains, we derive To oc 2 dN Ii. 

The packing fraction <pp where percolation of mobility 
domains occurs is lip = N N (1 — <frp/4>j) N . Mean- field 
percolation in configuration space occurs when the frac- 
tion of contacting mobility domains is l/(z — 1), where z 
is the average number of mobility domains that contact 
any single one in the <fi — > limit [3]. Since at percola- 
tion the average number of mobility domains that contact 
any single domain is 2 dN Hp, then 2 dN Hp/z = l/(z— 1). 
Since z increases with N, we predict (4>j — <j>p) oc TV -1 for 
large N. Thus, 4>p < <f>j for finite N whereas <\>p = <f>j as 
N — > oo. The dependence of <f>p on system size N sug- 
gests a diverging length scale £ in real space. For a system 
of hard spheres at packing fraction <j>, relaxation can only 
occur in subsystems of size Af ~ £ d where <j>p(N) > <f>. 
Using (4>j — <ft) oc 1/M we predict I oc (cf>j — <j))~ 1 / d . 

We can now express r Q in terms of <p. Since tq oc 
2 -dJV n -i K exp[A0j/Oj-0)] for large iV [H, the pack- 
ing fraction dependence of r Q can be determined from 
Eq. ©. For q£ < 1, 



q r a oc 



exp 



exp: 



A 4>.j 
<t>j-4 



for 
for 



- <PP, 
(3) 

where A and B are positive constants. In the large- 
TV limit when (j)p = 4>j the model predicts a Vogel- 
Fulcher divergence at <f>j. In this limit, the functional 
form and location of the divergence have been verified in 
experiments of hard spheres ]jj 21 1 • For finite N when 
(f)p < 4>j there is a power-law divergence for </> near <f>p 
and Vogel-Fulcher behavior far from cj>p. 

Finally, we emphasize that there are key differences 
between disordered hard-sphere systems in- and out-of- 
equilibrium. The largest connected cluster of mobil- 
ity domains has greater entropy than all other clusters, 
and equilibrium relaxation is described by the arguments 
above. However for large quench rates, the system is 
confined to a network of mobility domains that is de- 
tached from the percolating cluster. This results in a 
non-equilibrium glass state with <j> < (j>p, where is set 
by the size of non-percolating clusters, and the location 
of the glass transition depends on system preparation. 

B. Finite energy barriers: In contrast to hard 
spheres, activation is important in systems with finite 
energy barriers. We now extend the percolation model 



to include activated processes in systems at constant tem- 
perature T (with ks — 1). We again assume that only 
disordered states exist, as in frustrated geometries such 
as the pyrochlore lattice, polydisperse colloidal suspen- 
sions, and metallic glasses above the critical quench rate. 

For systems with finite energy barriers, the transition 
from glass to liquid is described by the percolation of 
bonds between local energy minima. Configuration space 
can be decomposed into basins of attraction surrounding 
each local minimum, and every point in configuration 
space can be mapped uniquely to a single basin 0]. At 
short times the system is confined to a basin, whereas at 
long times it will hop from one basin to another. Com- 
plete structural relaxation occurs once the system's tra- 
jectory percolates configuration space. 

To calculate r a , it is necessary to specify all possible 
percolating bond networks on which the system can relax 
and then select the subset that minimizes r a . To build 
the ensemble, we prescribe a maximum energy barrier 
height nT and draw bonds between minima with barriers 
below nT. This procedure defines networks of connected 
energy minima and if a network percolates it is included 
in the ensemble. For a given n, the fraction of bonds b(n) 
and the average time tq(u) to make transitions between 
two basins are given by 

b(n)= / P b {E)dE, (4) 
Jo 

T (n) oc bin)- 1 / P b (E)exp(E/T)dE, (5) 
Jo 

where Pb{E) is the distribution of energy barriers. For 
sufficiently large N, we can describe the properties of 
the percolating networks using a mean-field description 
and the only relevant variable is b(n). For a given n, a 
percolating network exists if b(n) is larger than a critical 
value bp, and the alpha relaxation time is 



T a (n) oc 



T (n)q 2 [b(n) - b P ] 
T (n)q- 6 



for q£ < 1, 
for q£ > 1. 



(0) 



Near percolation, when relaxation becomes slow, a 
saddle-point approximation holds and the system selects 
the n* that minimizes r a (n). Minimizing Eq. ([6]) for 
ql; < 1 gives 



b(n*) - bi 
b(n*) 



Ce n '/T (n*) 



(7) 



where C is the proportionality constant from Eq. (JSJ) with 
units of time. Eq. (JT]) can be solved to determine n* [22j]. 
Since Ce n " > r (n*) for all T > 0, Eq. (0 predicts that 
b(n*) > bp and there is no thermodynamic glass tran- 
sition at T > for systems with finite energy barriers. 
Furthermore, the T-dependence of r a (n*) depends on the 
distribution of energy barriers. 
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As in the case of hard spheres, our theory predicts 
that structural correlations such as the ISF or spin- 
spin correlations exhibit stretched-exponential relaxation 
with 1/3 < (3 < 1. These limits for the exponent (3 
agree with experimental values for many glass-forming 
materials From our percolation model, (3 increases 
slowly with time, but if it is measured near r a we ex- 
pect (3 = 1/3 for q£ » 1 and (3 = 1 for q£ <C 1, where 
£ oc [b(n*) — 6p]~ 1//2 . Our results for the limiting values 
of (3 are consistent with measurements at different £ val- 
ues in Lennard- Jones [23J] and magnetic 2J, |25j glasses 
and different q- values in magnetic [26[ and molecular 27 1 
glasses. 

The techniques in this section can be extended to treat 
hard spheres beyond the mean-field approximation. In- 
stead of energy barriers, hard spheres possess entropic 
barriers arising from the reduction of the effective di- 
mension of mobility domains far from CJ states. For a 
system that on average explores mobility domains at a 
distance R from the CJ state before making a transition 
to a new mobility domain, there is an average transi- 
tion time To (R) and a corresponding percolating network 
characterized by bond fraction b(R). In the large- N limit, 
the system will choose the effective domain size R* that 
minimizes the relaxation time. 

We have introduced a percolation model that gives 
rise to stretched-exponential relaxation in glass-forming 
materials. Stretched-exponential relaxation is commonly 
understood by assuming an underlying heterogeneity, ei- 
ther using trap models with distributions of waiting times 
and jump sizes [28| or dynamical heterogeneity in real 
space [291 ] . Percolation in configuration space naturally 
provides an origin to these heterogeneities via the struc- 
ture of the fractal percolating network, which is organized 
into densely connected blobs below £ and homogeneous 
nodes above £ |30j . Studies of the geometric properties of 
mean-field percolating networks, and dynamics on these 
networks, will lead to new predictions for slow dynamics 
and cooperative motion in glass-forming materials. 
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